WU-HEP-02-3 
TWC-02-1 

Comparison of space-time evolutions of hot/dense matter in v /JJvt7=17 and 130 
GeV relativistic heavy ion collisions based on a hydrodynamical model 

Kenji Morita, 1 '!] Shin Muroya, 2 '[] Chiho Nonaka, 3 '| and Tetsufumi Hirano 4, | 

1 Department of Physics, Waseda University, Tokyo 169-8555, Japan 
2 Tokuyama Women's College, Tokuyama, Yamaguchi 745-8511, Japan 
3 The Institute of Physical and Chemical Research (RIKEN), Wako, Saitama, 351-0198, Japan 
> ^Physics Department, University of Tokyo, Tokyo 113-0033, Japan 

o 

Based on a hydrodynamical model, we compare 130 GeY/A Au+Au collisions at RHIC 

(N ■ 

and 17 GeV /A Pb+Pb collisions at SPS. The model well reproduces the single-particle 

Oh" 

D ' distributions of both RHIC and SPS. The numerical solution indicates that huge amount of 

go ■ 

04 , to make a high temperature matter; longitudinal extent of the hot fluid in RHIC is much 



O 
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C^) • larger than that of SPS and initial energy density of the fluid is only 5% higher than the one 

in SPS. The solution well describes the HBT radii at SPS energy but shows some deviations 

o 
in , 
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from the ones at RHIC. 



INTRODUCTION 



One of the main purposes of ultra-relativistic heavy ion collision experiments is to explore the 



^ property of the hot and dense matter ]!]] . Recently the new experiment has begun to work at BNL- 

RHIC of which higher collision energy than other experiments up to now provides us chances to 
produce a new state of matter, quark-gluon plasma (QGP) with distinct possibility. However, the 
complicated collision processes composed of multiparticle productions and many-body interactions 
make it difficult to understand the properties of the hot matter. Therefore, a simple dynamical 
description of the system as a basis for the deeper understanding is indispensable. 

Relativistic hydrodynamical models are well-established phenomenological tools for describing 
high energy nucleus- nucleus collisions and subsequent multiparticle production |2|, |j| || ||, 0, || . 
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In this paper, we use a (3+l)-dimensional hydro dynamical model || with a QCD phase transition. 
We assume a cylindrical symmetry to the collision dynamics. Thus, our discussion is limited to the 
central collisions only. By virtue of the simple picture of our model, we can easily analyze both SPS 
and RHIC data with the same numerical code. Most of hydrodynamical calculations for RHIC data 
use Bjorken's scaling solution for the longitudinal direction. For example, Kolb et al. analyzed 
hadronic tranverse mass spectra and anisotropic flow 110]. Zschiesche et al. [11] investigated the 



HBT radii. These calculations assume the longitudinal boost-invariant infinite source. Though 



recently some hybrid models are used [[j^, [L^, 14 1 for the description of the hadronic phase, we here 



use a conventional description in which the hadronic phase is in local equilibrium. 

In this paper, concentrating our discussion on the central collisions, we reproduce the single- 
partcle spectra of hadrons at the beginning. In the hydrodynamical model, single-particle distri- 
butions are used as inputs in order to determine initial parameters rather than outputs. However, 
it is not trivial whether we succeed to reproduce experiments with "natural parameters" or not. 
Based on the solutions of hydrodynamical equations, we discuss the physical properties and the 
space-time evolution of the fluids in SPS and RHIC. We also evaluate the two-pion correlation 
functions and analyze the HBT radii. As a subsequent work of Ref. IE], we focus our discussion 



on comparison of the RHIC results and the SPS results based on the same numerical code. 

As is well-known, the two-particle correlation function gives us information on the size of the 



particle source [16, 17|. In the cases of the relativisitic heavy ion collisions, the correlation function 
tells us about the freeze-out which should be far from the static source. Thus, dynamical models 
such as hydrodynamical models are indispensable for understanding the relation between observed 
correlation functions and the space-time history of the system. However, up to now, any dynamical 
model assuming QGP failed to explain the experimental HBT radii in RHIC consistently with the 



single-particle spectra JllJ, |18|] and elliptic flow []19[, as being known as "HBT puzzle". We study 



the HBT radii in the framework of a hydrodynamical model which takes account of both transverse 
and longitudinal flow appropriately with a simple initial condition. 



In the next section, we explain our model. In Sec. Ill, we discuss the space-time evolution of 



the fluid. In Sec. IV, we present the result of two-particle correlation. Section ^ is devoted for 



the concluding remarks. 



1 In Ref. pel , a hadronic rescattering model is shown to reproduce these quantities. 
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II. HYDRO DYNAMICAL MODEL 

Let the system achieve the local thermal and chemical equilibrium shortly after a collision of 
two incident nuclei. This relaxation process cannot be described by the hydrodynamical model. 
The hydrodynamical model starts at initial time, ro, at which thermal and chemical equilibrium 
are established at least locally. The hydrodynamical equations are given as 

d^ v = 0. (1) 

We assume the perfect fluid for simplicity. Hence, energy-momentum tensor is given as 

= (e + P)U^U U - Pg^, (2) 

with U^, e and P being four velocities of a fluid element, energy density and pressure, respectively. 
These are treated as local quantities. We numerically solve the above equations together with the 
net baryon number conservation law, 

d^(n B Un = 0, (3) 

where n-Q is the net baryon number density and is also treated as a local quantity. Putting the z 
axis as a collision axis, we use a cylindrical coordinate system as follows; 

t = t cosh 7], (4) 

x = rcoscj), (5) 

y = r sin 0, (6) 

z = t sinh?7. (7) 

Focusing our discussion on central collisions, we may assume the cylindrical symmetry on the 
system. Therefore, by virtue of an identity U^U^ = 1, the four velocity can be expressed by two 
rapidity- like varibles Yl and It; 

U T = cosh(Y L - V) cosh F T , (8) 
U r > = sinh(Y L - 7]) cosh Y T , (9) 
U r ' = smhY T . (10) 

Most of hydrodynamical calculations which analyze RHIC data use Bjorken's scaling solution 
Yl = V- Putting the solution as an ansatz reduces numerical tasks very much but the analyses are 
limited to midrapidity region only. We solve not only transverse expansion but also the longitudinal 
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expansion explicitly. Numerical procedure for solving the coupled equations (||) and (||) is explained 
in Q. In this algorhythm, we solve the entropy and baryon number conservation law explicitly. 
Throughout our calculation, the total energy, entropy and baryon number are conserved within 5% 
of accuracy at the time step 5r = 0.01 fm/c. 

In order to solve the hydro dynamical equations, we must fix the equation of state (EoS). We 
adopt a bag model EoS with a first order phase transition. The QGP phase is composed of a free 
gas of massless u, d, s quarks and gluons. Hadronic phase is also assumed to be a free gas but with 
excluded volume correction. All hadrons are included up to 2 GeV/c 2 of mass except for hyperons. 
Putting the critical temperature as T c = 160 MeV at vanishing baryon density, we get the bag 
constant B 1 ^ = 233 MeV. We display the pressure as a function of temperature and baryonic 
chemical potential in Fig. |]. See [^] for the further detail of the EoS and numerical treatment of 
the first order phase transition in solving the hydrodynamical equations. 

We assume that the system achieves local equilibrium and begins to expand hydrodynamically 
at t = To = 1.0 fm/c. We put the initial conditions on this hyperbola. Bjorken's scaling solution 
is used as the initial condition of the longitudinal flow. Transverse flow is simply neglected at 
the initial time. We parameterize the initial energy density distribution e(ro,r],r) and net baryon 
number density distribution n^(ro,r], r) as, 



e(r ,r/,r) = e max exp 
nB(ro,V, r ) = n B0 \ exp 



2 — 0(M - Vo) 5 — 0(r - r ) 



2-^ 

(jl - VdY z 
2-al 



+ exp 



2-a\ 
2-al 



exp 



ro) 



2-al 



(r - r Q ) 



(11) 
(12) 



The energy density distribution of the longitudinal direction ( [TT[ ) has a central plateau characterized 
by r/o and a gaussian tail whose width is given by (Fig. ^), while the net baryon number 
distribution is a superposition of the two gaussians of which peaks exist at ±770 ■ For the transverse 
direction, both are parametrized by a flat region with gaussian smearing near the edge (Fig. For 
a nucleus with mass number A, the relation among these quantities is given by a r + tq = f.2^4 1 / 3 . 2 
Once these parameters are fixed, we can solve the hydrodynamical equations and pursue the space- 
time evolution of the fluid. These initial parameters are so chosen that the model reproduces the 
single-particle spectra measured in the experiments. The single-particle spectra can be calculated 



by making use of of the Cooper- Frye formula [22| 

dNi gi 



d 3 k (2tt) 



1 



K » a ° exp[(C/^- / , B )/T f ]Tl' 



(13) 



2 We adopt the initial condition as a natural and the simplest extension of the (1+1) -dimensional Bjorken's picture 
M and as a basis for the further improvement. 
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where gi is a degeneracy of the hadrons and Tf is a freeze-out temperature. The sign is plus for 
fermions and minus for bosons. Integration is performed on 3-dimensional freeze-out hypersurface 
S. By virtue of the Lagrangian hydrodynamics, contribution from the time-like hypersurface is 
expected to be small and the space-like hypersurface dominates the particle emission at freeze-out; 
we employ the non-covariant prescription k^da^ ~ k T da T for the sake of simplicity in the numerical 
treatment. In this approximation, total counted energy evaluated from Eq.(|l3|) is slightly larger 
than 90% of the total energy of the initial fluid; we regard the approximation works well enough. 
At SPS, we assume that the freeze-out occurs at a energy density ef and at a temperature Tf, at the 
RHIC energy. We also assume that the thermal and the chemical freeze-out are taken to happen 
simultaneouly. We show the freeze-out lines and the phase boundary on T — /xb plane in Fig. ||. 
Note that two freeze-out lines in the figure do not differ at low baryonic chemical potential (Fig. Q). 

We take account of the particles emitted from resonance decay as well as the direct emission 
from the freeze-out hypersurface. We include the decay processes p — ► 2tt, uj — ► 3n, rj — ► 3tt, 



K* — » ttK, and A — * Nir [g^, 24]. These resonances are also assumed to be thermally emitted 
from the freeze-out hypersurface. 

Two sets of initial parameter are summarized in Table |. Figures [B||7] show single-particle 
spectra in 17 ^4GeV Pb+Pb collisions at SPS. Our model well reproduces the experimental data 
with parameters in Table |[ Also in 130 AGeV Au+Au collisions, our model shows good agreement 
with the data as in Figs. B-lllj. However, we note that our model fails to produce enough number of 
anti-protons and overestimates the kaon yield in Fig. [H], where we multiply factors of 0.6 for kaons 



and 3.5 for anti-protons for clear comparison of the slopes [15]. This discrepancy may indicate the 
need for more sophisticated freeze-out mechanism. 

III. SPACE-TIME EVOLUTION 

In this section, we present the numerical solution of the hydro dynamical equations and discuss 



differences in the space-time evolutions of the fluids between RHIC and SPS. Figures 12 and 12 show 



space-time evolution of the fluid on the transverse plane. We also display the space-time evolution 



on the longitudinal plane in Fig. 14. Longitudinal flow and transverse flow are shown in Figs. 15 
and [l6|, respectively. From Table |, the maximum energy density in RHIC is only 5% higher than 
the one in SPS. Only 5% higher energy density for the almost 50% larger dN/dY seems surprising 
result. We show the number density of the thermal negative pions emitted into midrapidty as 
a function of the space-time rapidity rj of the freeze-out point (Fig. [1?]). This figure informs us 
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that thermal contribution of the volume element at r\ = to the particles into midrapidity is 
only 9% larger in RHIC than in SPS. However, the wider region of the freeze-out hypersurface 
in 7] contributes to the midrapidity particle distribution in RHIC more than in SPS. As a result, 
factor 1.5 times larger number of particles obtained in the midrapidity region after summing up 



particles emitted at different rj. The difference between RHIC and SPS in Fig. 17 originates in the 
longitudinal extent rj^ + a^ (see also Fig. |2|) and is direct consequence of longitudinal dynamics. We 



also plot the entropy per unit flow rapidity dS/dY[^ in Fig. 18. This is a conserved quantity if the 
boost-invariance is kept. In both RHIC and SPS, reflecting the deviation from the scaling solution 
shown in Fig. entropy is shifted to the larger flow rapidity. Reduction of entropy at Yl = 0, 
where = r] = always holds, comes from dY^/dr] which is larger than unity [25]. Thus, the 
shift at RHIC is smaller than the one at SPS since the deviation from the boost-invariant solution 
is small (Fig. |l5|). At SPS, the difference of dS/dY^ between the initial stage and the final stage is 
larger than the case at RHIC. 

Though our maximum energy density 6.0 GeV/fm 3 at initial is also so smaller than other 



calculation |JXC|, |llj, 26], this is due to the difference of initial time and transverse energy density 
profile of which nuclear thickness is considered. As for the average energy density at midrapidity, 
we get (eRHic) = 3.9 GeV/fm 3 and (esps) = 3.77 GeV/fm 3 . (epjuc) is a little smaller than an 



estimation of Ref. [27|. As a result of such a little difference in energy density, the space-time 
evolutions of two cases do not alter much in Figs. ^ and The most different point is a 
longitudinal extension of the fluid, t]q + a v . In RHIC, it is twice as large as in SPS. This is a 
consequence of much higher collision energy at RHIC. Indeed, the total energy of the fluid is 25290 
GeV at RHIC, which is 99% of total collision energy. Hence, higher collision energy does not lead 
to higher energy density but is used to produce the matter with large volume at To =1 fm/c. 

The output from the fluids is summarized in Table ||. Total net baryon number of the fluid 
is much smaller in RHIC than in SPS, as well as the mean chemical potential on the freeze-out 
hypersurface. This difference can be seen in the space-time evolution of temperature onij-r plane 
(Fig. |i~4| ). As shown in Fig. ||, phase boundary can no longer be specified by temperature only 
but depends on both temperature and chemical potential in high net baryon density. For example, 
T = 158 MeV corresponds to the hadronic phase at vanishing baryon density. However, it can be 
in QGP phase at /^b = 400 MeV and be in mixed phase at some /^b- This behavior is seen T = 158 



MeV contour in Fig. 14, where the baryonic chemical potential becomes higher as r] increases and 
a fluid element near r\ ~ 0.5 stays at mixed phase for a long time. Such behavior does not appear 
in the RHIC case where chemical potential is small enough to characterize the mixed phase by the 
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almost constant temperature. 



In Fig. 15 where deviation from Bjorken's scaling solution Y^ — r] is plotted, acceleration is 
larger in SPS than in RHIC because of the steeper pressure gradient of 77-direction. Finally, the 
lifetime of each phase is also shown in Table ||. 



IV. TWO-PARTICLE CORRELATION 

In this section, we present the result of the two-pion correlation function and HBT radii based 
on the numerical solution of the relativistic hydrodynamical equation. For simplicity, we assume 
that all the pions are emitted from a chaotic source and neglect the resonance contribution. Then, 
the two-particle correlation function is easily calculated through 

where = {k^ + k^)/^, q^ = k^ — k%, respectively p8| , [29[| . Here is on-shell mometum of i-th 
pion. We put 

I(q fi ,K^)= J K T da T (xWf(h,x)f(k 2 ,x) e^ x " , (15) 

so that 1(0, k^) reduces to the Cooper-Frye formula with f(k,x) being the Bose-Einstein distribu- 
tion function. Considering the experimental momentum acceptance, we integrate the correlation 
function with respect to the average momentum in region 0, as 

CW\ -1 1 fnKrdK T dY\I( q ^,Kn\ 2 

° 2{q )|n " 1 + f n K T dK T dYI(0,k»)I(0,k»y (16) 
The HBT radii are obtained by fitting the calculated correlation function (|I7J) to gaussian fitting 
function 

C2fit(<? M ) = 1 + exp(-itlg ide <? side - -Rout^out - ^longftong _ ^ol<?out<7long) • (17) 



For RHIC data, in which rapidity acceptance \Y\ < 0.5, the out-long cross term R Q \ [3C] can be 
ignored. According to the azimuthal symmetry, we can put Ki = K x , q s id e = % and (/out = Qx- 



Results of Kf dependence of the HBT radii are presented in Fig. [19j and |20|, where we show the 
transverse momentum dependence HBT radii of the SPS Pb+Pb collisions and Mp = \J + m% 
dependence of the RHIC Au+Au collisions, respectively. In addition to the three radius parameters, 
we also present the ratio of i? out to i? s id e for better comparison between two collisions [plf] . 

Sideward HBT radii (upper figures in Figs, [ll] and 20) are consistent with the experiments in 



both RHIC and SPS. Larger radii than other calculations come from initial large transverse size of 
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the fluid. Outward HBT radii show quantitative agreement for the SPS data. However, qualitative 
behavior shows some deviation from the experimental data; our result takes the maximum value at 
Kt ~ 0.3 GeV/c while the experiment data seem to monotonically decrease except for the highest 
Kt bin. For the RHIC data, experimental data show steep decrease with Mj. On the other hand, 
our results are similar to the ones of SPS because of the similarity in the space-time evolution of 
both fluids. As for the longitudinal HBT radii, our model reproduces the qualitative behavior of 
the results of both experiments but shows a little overestimate at low Mt of RHIC result. Our 
result suggests that the longitudinal finite size effect is essential for understanding the behavior of 
-Riong even at RHIC because other calculations assuming infinite boost invariant source show larger 
i?i ong [11, 19]. HBT radius in the longitudinal direction has two kinds of origin; spatial extent of 



the fluid [30 1 and thermal suppression caused by rapid expansion. The emission region is roughly 
characterized by a product of the Boltzmann factor exp[— wit cos1i(Yl — Y)/T] and a geometrical 
factor (e.g., exp(— r] 2 /2(Ar]) 2 ) ). The deviation from the scaling solution causes stronger thermal 
suppression. As a result, our solution provides smaller longitudinal HBT radius. Our -R ut/-Rsidc 
moderately increases with in both SPS and RHIC. This tendency is also seen in a hybrid 
calculation |ll| in spite of the quite different description of the hadronic phase. Our results show 
very good agreement with the SPS result, while the RHIC data clearly shows different behavior. As 



concerns the behaviors of i? out and -R ut/-Rside) opaque source |32j is a possible explanation if the 
current formalism of the two-particle correlation is correct. Though the opaque property appears 
in the hydrodynamical model by virtue of transverse flow [j33j [, it is still insufficient to reproduce 
the small R ut m the RHIC experiment. More theoretical investigation will be required for solving 
the problem |4|. 



Figure 21 shows the correlation functions projected onto each components of relative momenta. 
Transverse momentum of emitted pair is 0.125 < < 0.225 GeV/c, which corresponds to the 



lowest momentum bin in Fig. 2C. In each figure, calculated correlation functions are corrected by 
a common factor as A = 0.6. The other kind of the reduction of the correlation function at qi = 
{i =side, out, long) is caused by an integration with respect to other components of the relative 
momenta. The resultant factor is in proportional to 1/Ri. Our correlation function seems to be 
consistent with the experimental data for the outward and the longitudinal direction in spite of 
the overestimation of the HBT radii. Because the HBT radii correspond to the inverse width of 
the correlation function, difference in width between R = 7 fm and R = 6 fm is only about 5 MeV 
in relative momentum. The small value of sideward correlation function at q s ide = indicates that 
outward and longitudinal HBT radii are larger than the experimental results. 
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V. CONCLUDING REMARKS 



In this paper, we investigate single-particle distributions and two-particle correlation functions 
in SPS Pb+Pb 17 AGeV collisions and RHIC Au+Au 130 ^GeV collisions based on a hydrody- 
namical model in which both longitudinal and transverse expansion are taken into account. As 
long as the single-particle spectra, the hydro dynamical model well describes both SPS and RHIC 
data. The initial parameter set in the model for both collisions indicates that initial energy den- 
sity in RHIC is only slightly higher than the one in SPS and much larger extent of hot matter 
is produced in RHIC, if we compare them at the same initial time and by similar initialization 
(Fig. H and Eqs. (11),(]12])). We have also discussed the space-time evolution of the fluids. Since 
the initial conditions are not different so much, temperature and transverse flow evolution do not 
show significant difference between SPS and RHIC. Only the equi-temperature contour of T = 158 
MeV shows the qualitative difference between SPS and RHIC due to the difference in net baryon 
number. Steeper pressure gradient in the longitudinal direction at SPS makes the deviation from 
the scaling solution larger than at RHIC. Two-pion correlation functions and the HBT radii are 
also investigated. Our model shows good agreement with the SPS data. For RHIC data, however, 
the outward and longitudinal HBT radii of our result are larger than the experimental data even 
though the dynamical longitudinal expansion and finite size effect are incorporated. 
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FIG. 1: Pressure P(T, /xb) distribution 
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FIG. 2: Initial energy density distribution at r = in redirection. Dashed line stands for the RHIC case 
while solid line stands for the SPS case. 
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for RHIC 130A GeV Au+Au 




FIG. 3: Initial energy density (left) and net baryon number density (right) distribution for RHIC and SPS. 
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FIG. 4: Phase boundary (thick line) and freeze-out lines on T— hb plane. Thin solid line denotes the constant 
energy density contour as a freeze-out condition in SPS. Dashed line stands for constant temperature line, 
which is a freeze-out condition in RHIC. 
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FIG. 5: Rapidity distribution of negatively charged hadrons in Pb+Pb collisions at SPS. Closed circles are 
experimental data which are taken from 1 35 . Solid and dashed lines stand for our result of total yield and 
contribution from direct particles. 
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FIG. 6: Rapidity distribution of net protons in Pb+Pb collisions at SPS. Meanings of symbols and lines are 
the same as in Fig. |^. 




FIG. 7: Transverse mass distributions of negatively charged hadrons and net protons in Pb+Pb collisions 
at SPS. Closed and open circles are the experimental data of negatively charged hadrons and net protons, 
respectively. Solid and dashed lines stand for charged hadrons and net protons of our result, respectively. 
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FIG. 8: Pseudorapidity rj p distribution of charged hadrons. Solid line shows our result (ir,K,p) including 
resonance contribution. Dotted line denotes contribution of the directly emitted particles from the freeze-out 
hypersurface. Closed circles are preliminary result from the PHOBOS Collaboration p6|. 
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FIG. 9: Transverse momentum spectrum of negatively charged hadrons. As in Fig. g, solid line and dot- 
ted line show total number of particles and directly emitted particles from the freeze-out hypersurface, 
respectively. Closed circles are data from the STAR Collaboration J37| . 
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FIG. 10: Transverse mass spectra of negatively charged hadrons. Solid line, dotted line and dashed line 
denote ir~ , K~ and p yield of our result. K~ and p spectra are scaled down by factor 0.1 and 0.01, 
respectively. Closed circles, open triangles and closed triangles are preliminary data from the PHENIX 
Collaboration 1081. 
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FIG. 11: Rapidity dependence of anti-proton to proton ratio. Experimental data are taken from the 
BRAHMS Collaboration |§. 
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FIG. 12: Temperature contour on r — r plane at SPS. Left: r\ = section. Right: 77 = 1.5 section. 
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FIG. 13: Temperature contour on r — r plane at RHIC. Left: i] = section. Right: 77 = 3.0 section. 
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SPS Pb+Pb, r=0.0 [fm] 



RHIC Au+Au, r=0.0 [fm] 
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FIG. 14: Temperature contour on rj — t plane at r = fm. Left figure shows the SPS case and right figure 
shows the RHIC case. Thick solid (dashed) line shows the phase boundary between mixed (QGP) phase 
and hadronic (mixed) phase. Freeze-out line is given as dotted line for the SPS. As for RHIC, Tf = 125 MeV 
contour corresponds to the freeze-out line. 
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FIG. 15: Space-time evolution ol longitudinal flow at r = 0. Deviation from Bjorken's scaling solution of 
longitudinal flow rapidity is plotted. Thick lines stand for the RHIC case and thin lines stand for the SPS 
case. In both cases, solid lines, dashed lines, dotted lines and dash-dotted lines denote r =3, 7, 11 and 15 
fm/c cases, respectively. 
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FIG. 16: Space-time evolution of transverse flow. Left and right column show the SPS and the RHIC cases, 
respectively. 
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FIG. 17: Number density of the particles emitted into the midrapidity region as a function of space-time 
rapidity of the source point. Solid line stands for the SPS case and dashed line stands for the RHIC case. 
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FIG. 18: Entropy per unit flow rapidity. In both figures (upper for SPS and lower for RHIC), solid lines 
stand for the quantities calculated on the freeze-out hypersurface and dashed lines stand for the one on the 
initial stage. 
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FIG. 19: HBT radii at SPS. From top to bottom, i? s ide, Rout, -Riong and R ut/Rside are plotted. Closed 
circles denote the experimental data from the NA49 collaboration [jjoj] . Open squares stand for our results. 
Experimental acceptance is 2.9 < Y < 3.4 in the laboratory system. 
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FIG. 20: HBT radii at RHIC. From top to bottom, i? s idc 7 -Rout, -Riong and -R ut/-fiWe are plotted. Closed 
circles denote the experimental data from the STAR collaboration |)ti| . Open squares stand for our results. 
Experimental rapidity acceptance is \Y\ < 0.5. Three data points correspond to 0.125 < Kt < 0.225 GeV/c, 
0.225 < K T < 0.325 GeV/c and 0.325 < K T < 0.45 GeV/c, respectively. 
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FIG. 21: Correlation functions projected onto each component of relative momenta. From top to bottom, 
sideward, outward and longitudinal correlation functions are displayed, respectively. Each correlation func- 
tions are integrated from to 35 MeV with respect to other two components and corrected by a common A 
factor. 
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Tables 



TABLE I: Initial parameter set. 



SPS Pb+Pb 



RHIC Au+Au 



Maximum initial energy 
density e max 

"Maximum" initial net 
baryon density neo 
Longitudinal gaussian 
width u n of initial energy 
density 

Longitudinal extension r/ 
of the flat region in the ini- 
tial energy density 
Longitudinal gaussian 
width <td of the initial net 
baryon density 
Space-time rapidity tju at 
maximum of the initial net 
baryon distribution 
Gaussian smearing pa- 
rameter a T of the trans- 
verse profile 
Freeze-out condition 



5.74 GeV/fm 3 
0.7 fm~ 3 
0.61 

0.48 

0.52 



0.82 



1.0 fm 



6.0 GeV/fm 3 
0.125 fm" 3 
1.47 

1.0 

1.4 



3.0 



1.0 fm 



E f = 70 MeV/fm 3 



T f =125 MeV 
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TABLE II: Output quantities from the numerical solutions. 



SPS Pb+Pb RHIC Au+Au 



Net baryon number 305 131 

Mean freeze-out tempera- 123.2 MeV 125.0 MeV 

ture 

Mean chemical potential 281.6 MeV 76.1 MeV 

at freeze-out (/xb) 

Mean transverse flow ve- 0.467c 0.509c 
locity (vt) of the fluid at 
hi < 0.1 

Lifetime of the QGP phase 2.67 fm/c 2.92 fm/c 

T QGP 

Lifetime of the mixed 11.12 fm/c 12.61 fm/c 
phase tmix 

Total lifetime of the fluid 17.74 fm/c 18.94 fm/c 
thad 



